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NORMAL BINARY HIERARCHICAL MODELS 


DANIEL IRVING BERNSTEIN AND SETH SULLIVANT 

Abstract. Each simplicial complex and integer vector yields a vector configuration 
whose combinatorial properties are important for the analysis of contingency tables. We 
study the normality of these vector configurations including a description of operations 
on simplicial complexes that preserve normality, constructions of families of minimally 
nonnormal complexes, and computations classifying all of the normal complexes on up 
to six vertices. We repeat this analysis for compressed vector configurations, classifying 
all of the compressed complexes on up to six vertices. 


1. Introduction 

Associated to a simplicial complex C with ground set [m] := {1,2, and an 

integer vector d G Z™ is an integral matrix Ac, a (defined in Section [2]). The hierarchical 
model associated to C, d is a log-linear model (i.e. discrete exponential family, i.e. toric 
variety) whose design matrix is the matrix Ac, a- The linear transformation represented 
by the matrix Ac,a takes an array u G computes lower order marginals of u 

according to the faces of the simplicial complex C. As for all log-linear models, important 
relevant problems are to study properties of the lattice ker^Ac.d, the polyhedral cone 
M>oAc,d := {Ac,ax : a; > 0} and the semigroup MAc.d := {Ac,ax : a; > 0,a; integral}. 
In the special case where d = 2 the vector of all twos, we call Ac, 2 a binary hierarchical 
model. 

A matrix A G is called normal if M>oA fl ZA = MA. For general matrices A, 

deciding membership in the set NA is an NP-complete problem. However, for normal A 
the problem is polynomial time solvable since membership in both M>oA and ZA have 
polynomial time solutions (via linear programming and Smith normal form, respectively). 
So developing methods to decide a priori whether or not a matrix is normal can be useful 
for various problems related to integer programming. In this paper we focus on the 
following fundamental normality question for the matrices of hierarchical models. 

Question 1.1. For what values of C, d is the matrix Ac,d normal? 

If Ac,d is normal and d' < d, then Ac,d' is also normal, so it makes sense to focus 
on the binary case of Question 11.11 as a first step. Part of our motivation for studying 
this problem comes from previous joint work of Ranh and the second author [^, where 
normality was identified as a key property of hierarchical models to be able to apply the 
toric fiber product construction to calculate a Markov basis. 

The special case where the underlying simplicial complex C is a graph and d = 2 was 
solved in 0 where it was shown that a binary graph model is normal if and only if the 
graph is free of K 4 , minors. A complete classification is also available for the special case 
where C is the boundary of a simplex and d is arbitrary j^. 

In this paper we prove some new results and collect known results that relate to the 
classification of normal hierarchical models. We also perform extensive computations to 
check normality of hierarchical models where C has a small number of vertices. Section [2] 
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gives a detailed introduction to hierarchical models, the construction and interpretation 
of the matrix ^.d, and reviews the dehnition of normality. Section |3] shows how the 
classihcation of unimodular complexes [ij extends to a classihcation of all normal binary 
models on m vertices whose simplicial complex has a facet with m — 1 vertices. Section H] 
reviews several operations on the simplicial complex that are known to preserve normality 
and also proves the new result that normality is preserved on taking a vertex link. Section 
|6] gives the same results, but for the property of being compressed. Section |5] explores 
constructions for minimally nonnormal complexes. Section [7] surveys our computational 
experiments to classify complexes that give normal and compressed hierarchical models 
on small numbers of vertices. 


2. Preliminaries on Hierarchical Models and Normality 

In this section, we explain how to construct the matrix ^c,d associated to a hierarchical 
model and cover background material on normality of affine semigroups. 

Definition 2.1. Let d = (di,..., d„) G Z"" be such that di >2 for each i. Then we dehne 
to be the vector space of all di x ■ • • x tables. That is, each u G is a table 

in n dimensions where the ith dimension has di levels. We denote entries in a table n 
by where (ii, ... ,in) ^ [di] x • • • x [d„]. We will use the shorthand i to denote 

[ill ■ ■ ■ 1 in)- 

Definition 2.2. Let C be a simplicial complex on ground set [n]. A facet of C is an 
inclusion-maximal subset F C [n] that is contained in C. We let facet (C) denote the 
collection of facets of C. 


Definition 2.3. Let d = (di,..., d„) G Z” and F = {/i, ..., f^} C [n]. We dehne d^? to 
be the restriction of d to the indices in F, i.e. d^ = {df^,... ,df^). 


Definition 2.4. Let C be a simplicial complex on ground set [n] and let d = (di,..., d„) G 
Z"' be such that dj > 2 for each i. Then we have the linear map 

7rc,d : ^ 0 

_FGfacet(C) 


dehned by 


u 




• ij ^ [dj]) '-t 0 


i^Efacet(C) 




Ui,,...,i„ : j G [d/J X • • • X [d/J 




/ 


F={fi,...,fk} V :iF=j 

The matrix Ac.d denotes the matrix in the standard basis that represents the linear 
transformation 7rc,d- 


We give an example of this construction. 


Example 2.5. Let C be the simplicial complex on ground set [3] = {1, 2, 3} with facet(C) = 
{{1},{2,3}}. The linear transformation 7rc,d maps a three-way tensor u = {uijk : i G 
[di], j G [d 2 ], k G [da]) G to the direct sum of a one-way tensor and a two-way tensor: 

Fc,di^) = ^ [di]) © CY'^^F ■ j ^ [d2], k G [da]). 

j,k i 
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Taking d = 2 = (2, 2, 2), the matrix Ac,d could be represented as follows: 


0 

1 

0 

1 

0 

1 

0 

u 

1 

u 

1 

u 

1 

1 

1 

0 

0 

0 

1 

0 

0 

0 

0 

1 

0 

0 

0 

1 

0 

0 

0 

0 

1 

0 

0 

0 

1 

0 

\0 

0 

0 

1 

0 

0 

0 

1/ 


the horizontal line dividing the image space into the parts corresponding to the two facets 
{1} and {2, 3} of C. 

It can be useful to be precise about how to write down the matrix Ac^d- The columns 
of the matrix Ac,d are indexed by elements of the set njG[n]['^j] rows are indexed 

by pairs {F,e) where F G facet(C) and e G The entry in position with row 

index {F,e) and column indexed by i G is 1 if e = {ij : j E F), otherwise it is 

equal to 0. We are interested in whether or not this matrix is normal. 

Definition 2.6. Let A G We use the following notation for the lattice, cone, and 

semigroup spanned by the columns of A: 

Zd := {dz : z G M>od := {dr : r G M>o} Nd := {dn ; n G M"}. 

Note that Nd C Zd n M>od. We say the h is a hole of d if h G Zd n M>od but h ^ Nd. 
If d is free of holes, we say that d is normal. 


Not much is known about normality of the matrices dc,d- To end this section, we 
survey some known results from the literature on normality of hierarchical models. 


Theorem 2.7. \ ln . Theorem 1] Let C be a graph and d 


and only if C is K 4 minor free. 


2. Then dc,d is normal if 


Theorem 2.8. Theorem 2.2] Let C he the simplicial complex whose facets are all 
m — 1 element subsets of [m]. Then dc,d is normal in precisely the following situations 
up to symmetry: 

(1) At most two of the dy are greater than 2. 

(2) m = 3 and d = (3, 3, d^) for any d^ G N. 

(3) m = 3 and d = (3,4,4), (3,4, 5), or (3, 5, 5). 

Theorem 2.9. Theorem f 7 and following remarks] LetC = [12] [23] [34] [14] be the four¬ 
cycle graph. Then dc,d is normal if d = (2, ^2, 2, d^) or d = (2, d 2 , 3, ^4) with ^2, d^ G N. 


3. Complexes with a Big Facet 

The classification of unimodular binary hierarchical models given in [l| easily extends 
to a classification of normality of all binary hierarchical models containing a big facet 
(defined below). Throughout this section we assume that d = (2,..., 2) and we abbreviate 
Ac,d by Ac. We also abuse language and say “C is normal/unimodular” to mean “dc is 
normal/unimodular.” 

Definition 3.1. A matrix d G Z*^^" is unimodular if the polyhedron PA,b = (x G : 
dx = 6, X > 0} has all integral vertices for every b G Zd fl M>od. 

Proposition 3.2. If A E is unimodular, then A is normal. 
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Definition 3.3. For any matrix A G we define the Lawrence lifting of A to be the 
matrix 

/A 0\ 

A(Al) = P ^1 ^ M(2*+bx2i 

where 0 denotes the s x t matrix of all zeroes and 1 denotes atxt identity matrix. The 
kernel of this matrix is 

{{u, —u) : u G ker(74)}. 

When a matrix M has a kernel of this form, we say that M is of Lawrence type. 

Definition 3.4. Let C be a simplicial complex on [n]. We define the Lawrence lifting of 
C to be the simplicial complex A(C) on [n + 1] that has the following set of facets: 

{[n]} U {F U {n + 1} : F is a facet of C}. 

In this case, we refer to the facet [n] as a big facet. Any simplicial complex F that contains 
a big facet F (i.e. a facet containing all bnt one vertex) is the Lawrence lifting of the 
complex link^(F) where v is the vertex of F not contained in F (see Definition 14.61 for 
the definition of hnk„(F)). In this case we say that F is of Lawrence type. 

The following proposition justifies the multiple definitions of “Lawrence lifting.” 

Proposition 3.5. Theorem 15] A{Ac) = AIa(c)- 

Proposition 3.6. Corollary 5] Let A be an integer matrix of Lawrence type. Then A 
is normal if and only if A is unimodular. 

Combining Propositions 13.51 and 13.61 gives the following theorem. 

Theorem 3.7. Let A(C) be a simplicial complex of Lawrence type. Then A(C) is normal 
if and only if C is unimodular. 

Proof. By Propositions 13.51 and 13.61 we know that A(C) must be unimodular. Then [I| 
Proposition 3.14 gives that A(C) is unimodular if and only if C is unimodular. □ 

See [I[ for a complete classification of unimodular C. Combining this with Theorem 13.71 
gives a complete classification of the normal simplicial complexes of Lawrence type. 

4. Normality Preserving Operations 

An essential ingredient in the classification results from [ll| and [I[ was a list of sim¬ 
plicial complex operations that preserve normality and unimodularity, respectively. This 
section describes the operations that are known to preserve normality. In particular, we 
list some results from Hi namely that normality is preserved by taking induced sub¬ 
complexes, contracting edges, and gluing two complexes along a face. We then provide 
a new normality-preserving operation - taking links of vertices. For completeness sake, 
we also describe adding cone vertices and ghost vertices which are easily seen to preserve 
normality. 

Definition 4.1. Let C be a simplicial complex on vertex set V and let L G C. We define 
the edge contraction to be the simplicial complex on {n} U D \ L by 

C/L := {S e C : S n L = (!>} U {{S \ L) U {v} : S n L 0}. 
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Definition 4.2. Let C be a simplicial complex on vertex set V and let v G V. We define 
the vertex deletion to be the simplicial complex on \ {n} by 

Lemma 4.3 ([HI, Lemma 2.1). Suppose that C is obtained from C by either 

(1) deleting a vertex (C = C\v) or 

(2) contracting an edge (C = CfL). 

Then if Ac,d is normal, so is Ae',d', where d' is obtained from d by either ( 1 ) deleting d^ 
where v is the deleted vertex or (2) deleting all di such that i E L and setting dy = minjgi di 
in the case of a edge contraction. 

Definition 4.4. Let C be a simplicial complex on vertex set V. If Ci,C 2 are simplicial 
complexes on subsets of V and S V, then we say that C has reducible decomposition 
(Cl, S', C2) if C = Cl UC2 and Ci nC2 = 2 ^. If C has a redncible decomposition, we say that 
C is reducible. 

Lemma 4.5 ([IH, Lemma 2.4). LetC be a reducible simplicial complex with decomposition 
(Cl, S, C2). Assume that di and d2 agree at all vertices in S. Let d be the vector obtained 
by identifying the corresponding entries of di, 6 . 2 . If Ac^,di, AC 2 A 2 normal, then .4,c,d 
is normal. 


Definition 4.6. Let S' G C be a face of C. Then the link of S' in C is the new simplicial 
complex 

hnk5(C) = {F\S :F eC aJidS^F}. 

When S = {n}, we simply write hnk„(C) ;= hnk{^}(C). 

Note that we can obtain link 5 (C) by repeatedly taking links with respect to vertices. 
That is if S' is a face of C and ffS > 2 and v E S, then 

links (C) = hnk^(links\{^}(C)). 

We show that taking links of vertices preserves normality after proving establishing a 
dehnition and proving two propositions. 


Definition 4.7. For v G we let denote the matrix that projects onto the hyper¬ 
plane orthogonal to v: 


Proposition 4.8. Let a E W be a nonzero vector. Let T be a sguare matrix with row 
11 ^ such that the remaining rows span the hyperplane {x : a^x = 0}. Then the following 
are true 


(1) T is invertible 

(2) Ta = Cfc 

(3) Pa = T-iPe,T. 


Proof. The hrst two properties are clear so we proceed to prove the third. Withont loss 
of generality, assnme k = 1 and so -jj^ is the top row of T. Since the hrst row |||jp of T 
is orthogonal to the rest, the hrst column of T~^ must be a. Note that the matrix Pe^T 
is T with the top row replaced by Os. Therefore 


T-lPe,T 


I - 


a ■ aT 
11 ^ 


2 
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which is the dehnition of Pa 


□ 


Proposition 4.9. Let A G 

is normal. 


he a normal matrix with columns Then Pa^A 


Proof. First we show that the proposition is true when A is of the following form 


( 1 ) 


/ p 

^n—l 




0 

\ 

B 

0 


where > 0. If = 0, then P^^ is the identity so Pa„^ = A which is normal. So assume 
> 0. Let (°) G XPa^A n M>oPa„^ (note that Pa„^ has zeros along the top row). So 

such that Bz = Br = b. We lift b to b G TLA fl M>ov4 
An-i)- Choose d G Z such that (t,z) — {t,r) + dtn > 0. 


there exist z G and r G 
as follows. Define t = (fi,.. 
Define 


h:=A 


so b G IjA. Define 


e : = 


(t, z) - (t, r) + dp 


and so b is in 


^>0 


A because 


h:=A 


Since A is normal, there exists n G N"" such that An = b. Then Pa„An = and so 
(°) is not a hole of Pa^A- So Pa„^ is normal. 

Now we let a„ be an arbitrary nonzero vector. If T is an invertible matrix, then 
TA is normal. To see this, note that the solution sets to the equations TAx = b and 
Ax = T“^b are identical. So b is a hole of TA if and only if T~^h is a hole of A. Let 
T be a matrix whose hrst row is |j ||2 and whose remaining rows are orthogonal to a^. 
Since T is invertible, TA is normal. Furthermore, Proposition 14.81 (2) implies that the 
last column of TA is e{. So by previous arguments, Pa„TA is also normal. Therefore, so 
is T~^Pa„TA. This is equal to Pa„A by Proposition 14.81 (3). □ 

Theorem 4.10. Assume Ac,d is normal. Let v be a vertex of C and let d' denote the 
vector obtained by deleting dy from d. Then Aiink„c,d' is normal. 

Proof. Let A = Aiink„c,d' and let 1 
remark, we can write 


•^c,d — 


Assume A G and B G a 

may assume that A, B, Ac,d all have full row rank. Let AT denote the matrix that results 
when we project Ac,d onto the subspace orthogonal to the last (d — l)n columns. By 



i;,d^ • 

By Lemma 2.2 in and the following 

'A 

0 

... O' 


0 

A 

... 0 


0 

0 

... A 


_B 

B 

... P 


so A( 


Without loss of generality, we 
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Proposition 14.91 A! is normal. Then if 1 < i < n, the ith column of A! can be expressed 
as the ith column of .4.c,d minus a linear combination of the last {d — l)n columns of .4.c,d, 
all of which are 0 in the top m rows. So this means that the top m rows of A! are 

{A 0 ... 0) . 

Furthermore, since A and v4.c,d both have full row rank, the hnal [d — l)n columns of 
Ac^d have rank {d — l)m + 1. Therefore A! has rank m. Since A also has rank m, we may 
delete the bottom {d — l)m + I rows of A! without affecting the rowspace and therefore 
normality. So the matrix (^A 0 ... O) , and therefore A, is normal. □ 

For completeness sake, we give two very simple normality-preserving constructions. 

Definition 4.11. Let C be a simplicial complex on vertex set [n]. Then we dehne cone(C) 
to be the simplicial complex on [n + 1] with the following facets 

{F U {n -I- 1} : F is a facet of C} 

and we refer to this complex as a cone over C with n -|- 1 as the cone vertex. 


Definition 4.12. Let C be a simplicial complex on ground set [n]. Let G{C) denote the 
same simplicial complex but on ground set [n -|- 1]. Note that the vertex n -|- 1 is not 
contained in any face of G{C). In this case we say that n -|- 1 is a ghost vertex. 


Proposition 4.13. Let C be a simplicial complex on [n] and let C he cone(C) or G{C). 
Let d he an integer vector of length n, let d' he an extension of d of length n -|- 1. Then 
Aic,d is normal if and only ifAc',d' is normal. 


Proof. If C = cone(C) then 


Aci,d' 


/ Ac,d 

0 



0 

Af-CA 


0 





V 0 


0 

/ 


If C = G{C) then 


Ac',d' — 


Both of these matrices are normal if Ac^d is. 


•^C,d) • 


□ 


As a hnal tool for proving normality of a complex, we describe a condition on the 
polyhedral cone M>o.4.c that allows us to remove a facet of a normal complex C and retain 
normality. To describe this construction, which is a special case of Lemma 2 in 
need an alternate representation of the matrix Aq. 


11 


we 


Definition 4.14 (Full-rank representation of Ac)- Let C be a simplicial complex on the 
set {1,2,..., 77,}. For each binary 77 -tuple b = 6162 ■ ■ - bn G (0,1}" we dehne the column 
vector Hb whose entries are indexed by the faces of C. For each face F of C, we dehne 
ab(F’) = 1 if restricting b to the indices given by F contains no Is; otherwise we dehne 
ab(F") = 0. The matrix A whose columns are the set {ab}be{o,i}" has the same toric ideal 
as Ac. In particular, A is normal if and only if Ac is normal. 

See 0 for more details about the above construction, including proof of correctness 
and generalization to arbitrary d. For the rest of this section, we let Ac denote the matrix 
described in Dehnition 14.141 We now illustrate this construction with an example. 
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Example 4.15. Let C = [1][23] on ground set {1, 2, 3}. Then the faces of C are {0, {1}, {2}, 


{3}, {2, 3}}. Then Ac can be represented 

as 

follows 



/I 

1 

1 

1 

1 

1 1 

1 \ 


1 

1 

1 

1 

0 

0 0 

0 


1 

1 

0 

0 

1 

1 0 

0 


1 

0 

1 

0 

1 

0 1 

0 



0 

0 

0 

1 

0 0 

0 / 


Proposition 4.16. Let C be a simplicial complex and let S be a minimal non-face. Let 
V be the simplicial complex obtained from C by including S as a facet. We define B to be 
the matrix for the l-i-representation ofR^oAx); that is, B satisfies 

= {x G : Bx ^0}. 

There is exactly one column in B corresponding to the face S. If this column has entries 
contained in {0, +1, —1}, then ifV is normal, so is C. 

Proof. This is a special case of Lemma 2 in [Tl| . □ 

Example 4.17. For example, for the complex C with facets {!}, {2, 3}, the matrix B is 

/Of 0 0 0 \ 

1-10 0 0 
0 0 0 0 1 

^“ 0 0 0 1 -1 ■ 

0 0 1 0-1 

Vl 0 -1-1 1 / 

Applying Proposition 14. 161 we see that the complex C with facets {1}, {2}, and {3} is also 
normal. 


5. Minimally Non-Normal Binary Hierarchical Models 

Using normality- and unimodularity-preserving operations to obtain the classification 
results in nil and [I| required a complete list of the minimal non-normal and non- 
unimodular complexes. The list in [ll| was particularly nice; it included just one complex. 
The list in [l| included one infinite family and six other complexes. With our given list 
of operations (vertex deletion, edge contraction, reducibility, taking a cone) the list of 
minimally non-normal complexes appears to be even larger. This section describes our 
current knowledge of theoretical results towards describing this list. In particular, we re¬ 
port on two inhnite families of minimally non-normal complexes and a great many other 
complexes which may or may not fit into infinite families. 

If ^c,d is normal and d' < d, then Ac,d' is also normal. Therefore, a good hrst step 
towards a general classihcation of normal Acd would be a classihcation of the case where 
d = (2, 2 ,..., 2). We now restrict our attention to matrices Ac,d where d = (2, 2 ,..., 2) 
and so we abbreviate Me, (2,2,...,2) by Ac- Since the matrix Ac depends entirely on the 
simplicial complex parameter, we abuse notation and say “C is normal” when we mean 
“Ac is normal.” We give a list of all the complexes that we know to be minimally non¬ 
normal, a notion we make precise with the following definition. 

Definition 5.1. Let C be a simplicial complex. We say that a non-reducible complex 
C without cone and ghost vertices is minimally non-normal if C is not normal, but any 
edge-contraction, vertex-deletion, or link of C is normal. 
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By our general setup at the beginning of Section |2l the indexing set of the columns of a 
matrix Ac is [2] x [2] x • • • x [2] and the rows are indexed according to this set. However, 
since we are working in the binary case, it feels more natural to index with 0 — 1 strings 
than with 1 — 2 strings, so we instead index with the set {0,1}". This also allows us to 
use the notions of Hamming weight and Hamming distance in their natural form. 

In IfJ, the second author shows that the class of normal graphs is minor-closed with 
K 4 , the complete graph on four vertices, as the unique forbidden minor. Theorem 15.21 
below gives a family of non-normal simplicial complexes that generalize K 4 as a type of 
minimal forbidden minor. 


Theorem 5.2. Let C he the simplicial complex on the n + 1 vertices {0,1,... ,?7,} whose 
facets are the set {0,u} and all n — 2-dimensional sets of vertices that do not include both 
0 and n. Then Ac is minimally non-normal. 

Before proving Theorem 15.21 we must set some notation and prove two propositions. 

Definition 5.3. If a, 5 G are binary n-tuples, then we dehne the Hamming distance 
between a and b as follows 

dnifl^ b') .= 7 ^ bi^. 

That is, dnia, b) gives the number of indices in which a and b are different. We define the 
Hamming weight of a binary u-tuple a to be 

wnia) := #{i : a* = 1} = dH{a,0). 

That is wh{ci) is the number of nonzero entries in a. 

If i = (io, ii, ■ ■ ■, in) is a binary n + 1-tuple, define 

1 ''' . (^0; ■ ■ • iki • • • 1 in)i 

i.e. the binary n-tuple that results when we omit the fcth entry. If H is a collection of 
binary n-tuples, then we define 

Q\k 

Proposition 5.4. Given a collection B o/2"“^ binary n-tuples, the following two condi¬ 
tions are eguivalent: 

(1) For all a,b & B, d^ia, b) >2 

(2) The Hamming weight of all elements in B have the same parity. 

Proof. It is clear that ([2]) implies ([T]). So assume that B satishes ([T]). Then B'^^ consists of 
all 2"^“^ binary n — 1 tuples, as no two binary n tuples of the same Hamming weight parity 
can disagree only in the hrst position. Without loss of generality, assume (0,..., 0) G H. 
Now we show that all tuples in B have even hamming weight. For i G H, let denote 
the number of ones in the hnal n — 1 entries of i. If i G H has > 0, then for some j > 1, 
\j = 1. Let i' be the binary n tuple that agrees with i everywhere except i' = 0. Then 
i/\i _ pAi some i" G B. Since ^^(i, i') = 1, i" = 1 -|- \i{mod2). So either the hamming 
weights of i, i" are the same, or they differ by 2. By induction on k\, the Hamming weight 
of i" is even. Therefore, so is the Hamming weight of i. □ 

When A = Ac for some C, checking that some b G M>o^ is in ZH is easy. The following 
proposition shows that we only need to check that b has integer entries. 

Proposition 5.5. For any simplicial complex C, the matrix Ac has a full-rank submatrix 
with determinant 1. 
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gives 


Proof. We assume that Ac is presented as in Definition 14.141 Proposition 2.2. in 
that the rank of Ac is the number of faces in C. We now give a procedure for generating 
the desired submatrix. Let Fi ,..., be an ordering of the facets of C. Let 
denote the faces of that are not contained in any earlier facet. For each /j, choose the 
column of Ac whose index has zeros at exactly the vertices of /j. This submatrix is upper 
triangular with Is along the diagonal. □ 


We are now ready to prove Theorem 15.21 


Proof of Theorem \5.2[ We first prove the minimality claim. Deleting 0 or n yields dAn-i] 
deleting any other vertex yields A(A„_ 3 U Aq). The link about 0 orn is G(9A„_2); the link 
about any other vertex has reducible decomposition (clA„_ 2 , A„_ 3 , c}A„_ 2 ). Contracting 
an edge [On], [Ofc] or [kn] for 0 < fc < n yields 9A„_i; contracting any other edge yields 
A(A „_3 U Aq). All of these resulting complexes are normal. 

As stated in the beginning of this section, we can index the set of columns of Ac by 
the set of binary n + 1-tuples, {0, l}”'^^. We show that C is not normal by showing that 
the following vector is a hole 


b 



ie{o,i}"+i 


Recall that we can index each entry of Ui, b by a facet F oi C and a binary tuple whose 
entries are indexed by the vertices of F. The entry of Ui indexed by {F, i^) is 1 if the 
binary representation of i agrees with ii? at the indices indicated by F, and 0 otherwise. 
If F is an n — 2-dimensional facet (i.e. has n — 1 vertices), then for any fixed i^, there are 
exactly four indices i such that ai has a 1 at the entry corresponding to (F, i^?). So each 
(F, ij?) entry in b is 1, and in particular, an integer. For the edge [On], then for any fixed 
i[ 0 n], there are exactly 2"'“^ indices i such that ai has a 1 at the entry corresponding to 
([On], i[on]). So an entry in b corresponding to ([On],i[on]) is which is integral when 
n > 3. So when n > 3, b G M>oA fl and so b G M>oA fl TLA by Proposition 15.51 

We now show that b ^ NA. For the sake of contradiction, assume there exists some 
n G such that An = b. If ni > 0 then for each n — 2-dimensional facet F, ai adds 

Hi to exactly one entry indexed by a pair (F, ii?). Since each such entry is 1 in b, n must 
be a 0 — 1 vector. Furthermore, exactly entries of n must be 1 since there are exactly 
distinct {F,ip) for each fixed n — 2-dimensional facet F. So the vector n describes 
a collection of 2"^“^ binary n -f 1-tuples, which we denote F, i.e. 


Acn = y^aj. 

i&B 


Since C includes each F that is an n — 1-element subset of {0,..., n — 1} and each 
(F, ii?) entry of b is 1, we must have > 2 for every i,j G B. So Proposition 

15.41 implies that all the elements of BA have the same parity. So for i G F, we must 
have io = A + • • • + *n-i {mod2). Since C includes each F that is an n — 1-element 
subset of {l,...,n} and each (F, entry of b is 1, the same argument shows that 
in = A + ■•• + in-i {jnod2). So A = A for all i G F. But then the entries of b 
corresponding to ([On], (0,1)) and ([On], (1, 0)) must be 0. But we know that these entries 
are 2"'“^. □ 


We can use the forbidden-minor classification of unimodular simplicial complexes given 
in jl| to give some more examples of minimally non-normal complexes. Namely, if C is 
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normal, but minimally non-unimodular, then its Lawrence lifting A(C) is minimally non¬ 
normal. The minimally non-unimodular complexes include the inhnite family {dAn U 
{'i’}}n>i, and six other complexes. We know which of these are normal, and which are 
not. Computations in Normaliz show that the complexes P 4 , Ji, (dehned below) 
are normal. We can also see that dA^ U {n} is normal by Lemma 14.51 since dAn and 
{v} are normal (unimodular even - see 0). and dAn U {n} has reducible decomposition 

(aA„, 0 ,{n}). 

Proposition 5.6. Let C be among the following complexes. 

(1) Pi, the path on 4 vertices 

(2) Ji, the complex [12] [15] [234][345] 

(3) J* the complex [134] [235] [245] (this is the Alexander dual of Ji) 

(4) dAn U {n} for n>l. 

Then A(C) is minimally non-normal. 

Proof. Let u denote the vertex added when creating A(C) from C. If we delete u we are 
left with a simplex, which is normal. If we take the link of the u, then we are left with 
C which is normal. If we take the link or deletion of any other vertex w ^ u, we are 
left with a Lawrence lifting of a link or deletion of C, call it C. Since C is minimally 
non-unimodular, C is unimodular and so A(C') is normal. If we contract an edge that 
connects to u, then we are left with a simplex, with is normal. If we contract an edge that 
does not connect to u, then we are left with a Lawrence lifting of an edge contraction of 
C, call it C. We can check that C is unimodular and so A(C') is normal. □ 

The remaining three minimally non-unimodular complexes are not normal. The com¬ 
plex J 2 = [12] [235] [34] [145] is not normal but not minimally so - contracting the edge [14] 
gives A(c}Ai U {n}). The other two minimally non-unimodular complexes are minimally 
non-normal, as we show in the following proposition. 

Proposition 5.7. The following complexes are minimally non-normal: 

(1) [123] [134] [145] [125] [623] [634] [645] [625], the boundary of the octahedron 

( 2 ) [1234] [3456] [1256], the Alexander dual of the boundary of the octahedron. 

Proof. A computation in Normaliz [3| shows that the hrst complex is not normal. The 
matrix for the second complex is the matrix for the no 3-way interaction model with 
I"! = i ’2 = I’s = 4, which is not normal [^, Theorem 6.4]. All links and induced sub¬ 
complexes are unimodular [I| and therefore normal. Contracting an edge in the hrst 
complex gives either the cone over a square or the Alexander dual of A 2 U Ai, both of 
which are unimodular. Contracting an edge in the second complex gives a complex whose 
matrix is the matrix for the no 3-way interaction model with ri = r 2 = 4 and r^ = 2 
which is normal [^, Theorem 6.4]. □ 

6. Compressed Hierarchical Models 

A stronger property for a matrix A G is to be compressed. This condition guar¬ 
antees normality and can be easier to verify directly through polyhedral computations. 
As in the normal case, we show that taking induced subcomplexes, edge contractions, 
reducible decompositions, and links preserves the property of being compressed. 

To state the dehnition of compressed vector conhguration we hrst need the dehnition 
of the toric ideal associated to the conhguration A. See for more details on toric ideals 
and their Grobner bases. 
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Definition 6.1. Let A G and K[a;] := ]K[a;i,... ,x„]. The toric ideal associated to 
the matrix A is the ideal 


(x“ — x’' : M, n e M”, Au = Av). 

Definition 6.2. If A G is an integer matrix such that the corresponding toric ideal 
I A C ]K[xi,..., x„] is homogeneous. Then we say that A is compressed if the initial ideal 
in^{lA) is generated by squarefree monomials whenever -< is a reverse lexicographic term 
order. 


It can often be easier to use the following characterization of compressed matrices. 

Theorem 6.3 ([l^). Let A G Z^^"" be an integer matrix such that the corresponding 
toric ideal I a C ]K[xi, ... ,x„] is homogeneous and let B be the matrix of facet defining 
inegualities of M.^qA. In other words, B satisfies 

M>oA = {Ax : X > 0} = {i/ : By > 0}. 

Suppose that the entries in each column of B are relatively prime. Then A is compressed 
if and only if AB^ is a 0/1 matrix. 

As the following Proposition states, compressed-ness is a stronger property than nor¬ 
mality, and weaker than unimodularity. Hence, to classify the compressed models, we can 
restrict attention to normal models. 


Proposition 6.4. //A G is compressed, then A is normal. If A E Z'^^" is unimod- 
ular, then A is compressed. 


Proposition 6.5. The following operations on simplicial complexes preserve the property 
of being compressed of Ac,d-' 

• gluing two complexes along a common facet 

• passing to induced subcomplexes 

• adding or removing cone vertices 

• contracting edges. 


Proof. See 



□ 


We also add to the list of compressed-ness operations computing a link of a vertex. 

Proposition 6.6. Let A G Z*^^” be a matrix with columns ai,...,cLj. Assume that 
the entries of a„ do not have a common divisor. Let A' be the matrix consisting of the 
first n — 1 columns of Pa^A. Let (f : K[xi,..., x„] —)■ ]K[xi,..., x„_i] be the substitution 
homomorphism mapping x„ to 1. Then Ia' = A{Ia)- 

Proof. Note that u G kerz(A') iff Pa„A(“) = 0 which is true iff Au = for some k eTj 
which is true iff n — (0,..., 0, ky G kerz(A). Now let be a generator of Ia>- 

Since n — (0,..., 0, ky G kerz(A), we have x“^ — x^x““ G Ia without loss of generality. 
Then note x“^ — x““ = 0(x“^ — x^x““). □ 

Theorem 6.7. Assume Ac,d is compressed. Let v be a vertex of C and let d' denote the 
vector obtained by deleting the entry for v from d. Then A \\r^v ,,c d' is compressed. 
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Proof. Let A = ^unk^c.d' and let B = v4.c\w,d'- As in the proof of Theorem 14.101 we know 
that we can write 


'A 

0 ... 

O' 

0 

A ... 

0 

0 

0 ... 

A 


B ... 

B_ 


and that projecting onto the space orthogonal to the last {d — l)n colnmns and deleting 
unnecessary rows, then deleting columns of zeros, leaves the matrix A. 

Let ]K[x] be the ring containing I a, and K[x, y] be the ring containing Iac^- By 
Proposition 16.61 I a is obtained from Iac a by plugging in 1 for all the y variables. Since 
IA is the toric ideal for a hierarchical model, it is homogeneous. Consider any reverse 
lexicographic order -< on ]K[x]. Extend this to a reverse lexicographic order -<' on ]K[x, y] 
by putting the y variables ahead of the x variables in the variable order (not in the term 
order). This means that if x“^ — x“ € has initial term x“^ with respect to then 
for any v such that y’^ x“ — x“ y’^ G Iac a i then y’' x“ is the initial term with respect 
to -<'. 

Since Ac ,a is compressed, there is a Grobner basis element 
whose leading term y^^x“^ is a squarefree monomial and divides y^^x“^. Without loss of 
generality, we can assume that u~^ ^ 0. Indeed, any y'“^x“^ — x“ y^ G Iac a with vA = Q 
must also satisfy u~ = 0, and so we can reduce y’'^x“^ — x“ y^ with respect to such 
polynomials to obtain a polynomial where y'^^ is reduced with respect to the Grobner 
basis of I Ac a • 

Now apply the map 0 : ]K[x, y] —)• ]K[x] that sets all y variables to 1. Then the leading 
term of 0(y*^^x“^ — x“ y’^ ) = x“^ — x“ is x“^ which divides x“^. This implies that the 
initial terms in any Grobner basis element of Ia with respect to the reverse lexicographic 
order are squarefree. Hence, ^iink„c,d' is compressed. □ 

7. Computational Results 

In this section we summarize the results of our computations to classify the normality 
and compressed-ness of Ac for complexes with few vertices. In particular, we highlight 
the simplicial complexes whose normality and compressed-ness, or lack thereof, cannot be 
determined from our existing construction methods for normal complexes or our examples 
of minimally nonnormal complexes. In particular, we wrote a Mathematica code to gen¬ 
erate, up to symmetry, all simplicial complexes for which normality and compressed-ness 
cannot be ruled out by looking at induced subcomplexes, vertex links, edge contraction, 
cone vertices, or reducibility. From the resulting lists we excluded all simplicial complexes 
that are graphs (where normality is already classified by H) or that are of Lawrence type 
(where normality follows from the classification of unimodular complexes 0). Then we 
use Normaliz to decide normality of the remaining complexes in our lists. The results 
are summarized below. 

7.1. Four and Fewer Vertices. Every simplicial complex on 1, 2, or 3 vertices is normal 
and all of these follow from the construction methods described above. On m = 4 vertices, 
the normality or nonnormality of all simplicial complexes can be decided by the methods 
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listed above. Among all complexes on 4 vertices, the only two that are nonnormal are the 
following: 

(1) |12||13||14||23||24||34| 

(2) |123]|14]|24]|34] 

both of which are clearly minimally nonnormah The compressed complexes on four ver¬ 
tices are exactly the normal ones. 

7.2. Five Vertices. Applying the procedure described above to simplicial complexes on 
5 vertices produces 14 complexes which the techniques mentioned above do not decide 
the normality of. Of these 14 complexes, 6 were normal. These are 

(1) [12] [13] [245] [345] 

(2) [12] [134] [235] [345] 

(3) [123] [124] [135] [245] 

(4) [12] [134] [135] [234] [235] 

(5) [123] [124] [134] [235] [245] 

(6) [123] [124] [134] [235] [245] [345]. 

The hrst complex is normal because it is equivalent to the four-cycle with d = (2, 2, 2,4), 
so normality follows from Theorem 12.91 The last complex is normal because it is a 
unimodular complex, in particular it is the Alexander dual of Ai U A 2 . For the other 
complexes, normality was verihed using Normaliz. 

The remaining 8 out of 14 complexes are all nonnormal. These are 

(1) [123] [124] [135] [245] [345] 

(2) [123] [124] [125] [134] [234] [345] 

(3) [123] [124] [125] [134] [235] [345] 

(4) [12] [134] [135] [145] [234] [235] [245] 

(5) [123] [124] [125] [134] [135] [234] [245] 

(6) [123] [124] [125] [134] [135] [245] [345] 

(7) [12] [134] [135] [145] [234] [235] [245] [345] 

(8) [123] [124] [125] [134] [135] [234] [245] [345]. 

The second to last complex on this list is verihed to be non-normal as a special case 
of Theorem 15.21 The other complexes on the list were shown to be non-normal using 
Normaliz. On hve vertices, the minimal non-normal complexes consist of these eight 
simplicial complexes, plus the two complexes of Lawrence type 

(1) [1234] [125] [235] [345] 

(2) [1234] [15] [235] [245] [345]. 

We can use the results from 4-vertex complexes and Theorem 16.71 and Proposition 16.51 
to show that every normal complex on 5 vertices is compressed, aside from the 5-cycle 

[12] [23] [34] [45] [15], 

7.3. Six Vertices. We applied the above procedure to simplicial complexes on 6 vertices. 
This produced 80 complexes for which the techniques above do not determine normality. 
Using some other techniques, described below, we were able to determine normality and 
compressedness of all of these complexes. The results are described in the table below. 
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Complex 

Normal? 

Notes 

Compressed? 

El 

[1456] [1236] [2345] 

No 

0 

No 

1 

[2346] [16] [2345] [15] 

Yes 

0 

Yes 

m 

[2346] [16] [2345] [125] 

Yes 

AC 

Yes 

a 

[236] [16] [2345] [145] 

Yes 

C 

Yes 

m 

[2346] [16] [2345] [1235] 

Yes 

AC 

Yes 

B 

[3456] [126] [1345] [123] 

Yes 

AR 

Yes 

B 

[346] [126] [345] [125] 

Yes 

AR 

Yes 

B 

[1346] [126] [345] [125] 

Yes 

C 

Yes 

8 

[1346] [126] [1345] [235] 

Yes 

C 

Yes 

B 

[3456] [126] [1345] [1234] 

Yes 

c 

Yes 

B 

[1246] [1236] [2345] [1345] 

Yes 

AR 

Yes 

B 

[26] [16] [345] [245] [13] 

Yes 

AR 

No 

B 

[236] [16] [245] [235] [14] 

Yes 

AR 

No 

B 

[2456] [236] [16] [1245] [123] 

Yes 

AR 

Yes 

B 

[236] [16] [245] [235] [124] 

Yes 

AC 

No 

B 

[236] [16] [245] [145] [234] 

Yes 

C 

No 

IQ 

[2356] [2346] [16] [1235] [1234] 

Yes 

AC 

Yes 

B 

[346] [126] [345] [135] [123] 

Yes 

AR 

No 

B 

[2456] [136] [126] [2345] [123] 

Yes 

AR 

Yes 

B 

[3456] [2356] [1346] [126] [123] 

Yes 

AC 

Yes 

^1 

[136] [126] [345] [245] [124] 

Yes 

C 

Yes 

B 

[136] [126] [2345] [1345] [124] 

Yes 

AR 

Yes 


[1356] [346] [126] [1345] [124] 

Yes 

AR 

Yes 

BSl 

[2356] [1346] [126] [2345] [124] 

Yes 

C 

Yes 

Q| 

[3456] [1346] [126] [1345] [125] 

Yes 

AR 

Yes 


[2456] [1356] [1246] [1236] [125] 

Yes 

AR 

Yes 


[2356] [1346] [126] [2345] [1234] 

Yes 

AC 

Yes 


[3456] [1346] [126] [2345] [1234] 

Yes 

C 

Yes 


[3456] [1246] [1236] [1345] [1234] 

Yes 

AR 

Yes 

BSl 

[3456] [1246] [1236] [1245] [1235] 

No 

C 

No 

^1 

[26] [16] [245] [145] [23] [13] 

Yes 

C 

Yes 

B 

[456] [356] [246] [136] [126] [14] 

Yes 

AR 

No 


[246] [236] [16] [245] [235] [15] 

Yes 

AR 

Yes 


[236] [16] [235] [125] [234] [124] 

Yes 

AC 

Yes 

Q| 

[3456] [2456] [136] [126] [234] [123] 

Yes 

AR 

Yes 


[246] [136] [126] [245] [235] [123] 

Yes 

AR 

No 


[3456] [2456] [136] [126] [2345] [123] 

Yes 

AR 

Yes 


[2346] [1346] [126] [2345] [1345] [123] 

Yes 

AC 

Yes 


[2356] [1356] [2346] [1346] [126] [123] 

Yes 

AR 

Yes 

^1 

[136] [126] [235] [125] [134] [124] 

Yes 

AR 

Yes 

B!il 

[136] [126] [245] [235] [134] [124] 

Yes 

AR 

No 

B 

[136] [126] [345] [245] [134] [124] 

Yes 

C 

Yes 

IQ 

[3456] [2456] [136] [126] [134] [124] 

Yes 

C 

Yes 

IQ 

[136] [126] [145] [135] [234] [124] 

Yes 

AR 

No 

IQ 

[136] [126] [235] [135] [234] [124] 

Yes 

AR 

Yes 

IQ 

[136] [126] [345] [135] [234] [124] 

Yes 

AR 

No 

IQ 

[2346] [136] [126] [2345] [1345] [124] 

Yes 

C 

Yes 

IQ 

[3456] [2356] [1346] [126] [2345] [124] 

Yes 

c 

Yes 

IQ 

[2346] [1346] [126] [2345] [1345] [125] 

Yes 

c 

Yes 

IQ 

[246] [136] [126] [345] [245] [135] 

No 

c 

No 

^1 

[1356] [2346] [1346] [126] [1235] [1234] 

Yes 

AR 

Yes 

Continued on next page 
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Table 1 continued from previous page 



Complex 

Normal? 

Notes 

Compressed? 

51 

[3456] [2356] [1346] [126] [2345] [1234] 

Yes 

AC 

Yes 

52 

[2456] [1356] [1246] [1236] [1345] [1234] 

No 

C 

No 

53 

[2356] [1346] [1246] [1236] [2345] [1234] 

Yes 

AR 

Yes 

54 

[1346] [1246] [1236] [2345] [1245] [1235] 

Yes 

AR 

Yes 

55 

[2456] [1356] [1246] [1236] [2345] [1345] 

No 

C 

No 

56 

[346] [246] [136] [126] [245] [125] [14] 

Yes 

AR 

No 

57 

[456] [356] [256] [146] [136] [126] [15] 

Yes 

AR 

Yes 

58 

[356] [246] [136] [126] [235] [234] [123] 

Yes 

AR 

Yes 

59 

[356] [246] [136] [126] [345] [234] [123] 

Yes 

AR 

No 

60 

[346] [246] [136] [126] [345] [245] [123] 

Yes 

AR 

No 

61 

[136] [126] [135] [125] [234] [134] [124] 

Yes 

AR 

Yes 

62 

[246] [136] [126] [345] [245] [135] [125] 

No 

C 

No 

63 

[3456] [126] [2345] [1345] [1245] [1235] [1234] 

No 

C 

No 

64 

[2356] [1356] [2346] [1346] [126] [1235] [1234] 

Yes 

AR 

Yes 

65 

[3456] [2356] [1346] [126] [2345] [1345] [1245] 

No 

C 

No 

66 

[1456] [1356] [1246] [1236] [2345] [1345] [1234] 

No 

C 

No 

67 

[2456] [1356] [1246] [1236] [2345] [1345] [1234] 

No 

C 

No 

68 

[2456] [2356] [1346] [1246] [1236] [2345] [1234] 

Yes 

AR 

Yes 

69 

[456] [356] [246] [136] [126] [345] [234] [123] 

Yes 

AR 

No 

70 

[136] [126] [235] [135] [125] [234] [134] [124] 

Yes 

AR 

Yes 

71 

[346] [246] [136] [126] [345] [245] [135] [125] 

No 

C 

No 

72 

[1456] [1356] [2346] [1346] [126] [2345] [1345] [125] 

No 

C 

No 

73 

[3456] [2456] [2356] [1346] [1246] [1236] [2345] [1234] 

Yes 

u 

Yes 

74 

[2456] [2356] [1346] [1246] [1236] [1345] [1245] [1235] 

No 

c 

No 

75 

[3456] [2456] [2356] [1346] [1246] [1236] [2345] [15] [1234] 

No 

c 

No 

76 

[236] [136] [126] [235] [135] [125] [234] [134] [124] 

Yes 

AR 

Yes 

77 

[3456] [1456] [1356] [2346] [1346] [126] [2345] [1345] [125] 

No 

C 

No 

78 

[3456] [2456] [2356] [1346] [1246] [1236] [1345] [1245] [1235] 

Yes 

u 

Yes 

79 

[3456] [2456] [2356] [2346] [16] [2345] [1345] [1245] [1235] [1234] 

No 

0 

No 


U: Complex is unimodular. 

AC: Complex is proved normal using Proposition 14.161 and the fact that one face can 
be added to produce a cone. 

AR: Complex is proved normal using Proposition 14.161 and the fact that one face can 
be added to produce a reducible complex. 

C: Complex was proved normal, or not normal, by a computation in Normaliz. 

O: Normality/Nonnormality determined by other method or theoretical result. 

The complexes indicated by an O had their normality or nonnormality decided by a 
theoretical result mentioned in the paper. Complex 1 has the same matrix as Ac,d where 
C is the 4-cycle and d = (2, 2, 2, 8), which is normal by Theorem 12.91 Complexes 73 and 
78 are unimodular and therefore normal. Complex 0 has the same matrix as Ac,d where 
C is the 3-cycle and d = (4,4,4) which is not normal by Theorem 12.81 Complex 79 is 
an instance of the minimally non-normal family from Theorem 15.21 For complexes 55, 
65, 67, 72, 75, and 77 our computations with Normaliz were initially unsuccessful. The 
Normaliz team was able to conhrm nonnormality of these examples using the current 
development version. 
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Proposition 14.161 was used to prove normality in a number of cases where we added 
a face and either produced a reducible complex or a cone complex that was known to 
be normal. These are indicated in the table with an AR or AC respectively. With our 
expanded list of complexes known to be normal, we can apply this same technique over 
and over again until we no longer add any new complexes to the set that we know to be 
normal. This technique only gives us three more complexes - adding the facet [236] to 
complex 27 yields 37, adding [23] to 41 yields 56 and adding [134] to 20 yields 41. In all 
cases, the added face satisfies the conditions of Lemma 14.161 and so complexes 27,41, and 
20 are normal. 


8. Future Directions 

Results from the previous section suggest that a complete classihcation of complexes 
C such that Ac in normal is probably extremely complicated. At any rate, our current 
tools for proving normality or nonnormality are inadequate for completing a classification. 
This suggest the following problems: 

Problem 8.1. Develop new construction techniques for producing complexes C such that 
Ac is normal. 

The technique based on Proposition 14.161 remains mysterious and application of it in 
more settings (and generalizations of it) would require an understanding of the facet 
defining inequalities of M>o.4.c for families of simplicial complexes. 

Problem 8.2. Develop new methods for constructing holes in the semigroups N.Ac. 

In particular, it would be worthwhile to analyze the holes produced in the nonnormal 
examples found above and try to generalize them to find infinite families of minimal 
nonnormal complexes. 

While a complete classification of the complexes C such that Ac is normal is out of 
reach at present, it might be reasonable to try to classify low dimensional complexes that 
yield normal vector conhgurations. For example, the set of 1-dimensional C such that 
Ac is normal is extremely simple (the set of A' 4 -minor free graphs 0 ). so perhaps 2 - 
dimensional complexes C such that Ac are normal also have a simpler description. As a 
starting point, we propose the following: 

Problem 8.3. Classify the two-dimensional manifolds C such that Ac is normal. 

In principle, the solution to Problem 18.31 should yield a relatively simple list of com¬ 
plexes. Indeed, a 2-dimensional manifold with many triangles seems likely to have a minor 
isomorphic to K/^ or the complex [123] [14] [24] [34], so there seem to be few possibilities. 
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